Fast atomic transport without vibrational heating 
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We use the dynamical invariants associated with the Hamiltonian of an atom in a one dimen- 
sional moving trap to inverse engineer the trap motion and perform fast atomic transport without 
final vibrational heating. The atom is driven non-adiabatically through a shortcut to the result of 
adiabatic, slow trap motion. For harmonic potentials this only requires designing appropriate trap 
trajectories, whereas perfect transport in anharmonic traps may be achieved by applying an extra 
field to compensate the forces in the rest frame of the trap. The results can be extended to atom 
stopping or launching. The limitations due to geometrical constraints, energies and accelerations 
involved are analyzed, as well as the relation to previous approaches (based on classical trajecto- 
ries or "fast-forward" and "bang-bang" methods) which can be integrated in the invariant-based 
framework. 



PACS numbers: 03.65.Ta,03.65.Ca,03.65.Nk 



I. INTRODUCTION 

A key element to attain an exhaustive control of states 
and dynamics of cold atoms and ions is their efficient 
transport by moving the confining trap. In spite of the 
broad span of conditions, heating mechanisms, transport 
distances (from microns to tens of centimeters), trans- 
port times, and accelerations involved, there are some 
common elements and objectives that allow for a rather 
generic theoretical treatment as the one presented in this 
paper. Transport should ideally be lossless, fast, and lead 
to a final state as close as possible ( "faithful" ) to the ini- 
tial one, up to global phase factors, in the frame of the 
transporting trap. The later requirement is character- 
ized in the most demanding applications as a high-fidelity 
condition, or more generally as a no-heating or at least 
minimal heating condition; equivalently, by the absence 
or minimization of vibrational excitations at the end of 
the transport. Note that reaching a faithful final state is 
not incompatible with some transient excitation in the in- 
stantaneous basis at intermediate times, i.e., the process 
does not have to be slow (in the usual quantum mechani- 
cal jargon, "adiabatic"), although slowness is certainly a 
simple way to avoid heating, at least for ideal conditions. 

Efficient atom transport is a major goal for many appli- 
cations such as quantum information processing in mul- 
tiplexed trap arrays [l], Q or quantum registers [3| ; con- 
trolled translation from the production (cooling) cham- 
ber to interaction or manipulation zones [343; accurate 
control of interaction times and locations, e.g. in cavity 
QED experiments [7[ , quantum gates [8j or metrolog y [9 j ; 
and velocity control to launch [1(|, or stop atoms [llLll^. 

Different approaches have been implemented. Neutral 
atoms have been transported as thermal atomic clouds 
@, [HI , condensates [lj] , or individually [H, [l6| , using 
magnetic, or optical traps. The magnetic traps can be 



translated by moving the coils mechanically Q , by time- 
varying currents in a lithographic conductor pattern [llj , 
or on a conveyor belt with a chain of permanent magnets 
[l~8l ] . Optical traps can be used as optical tweezers whose 
focal point is translated by moving mechanically lenses 
4, 19], or by traveling lattices (conveyor belts) made with 
two counterpropagating beams slightly detuned one re- 
spect to each other [la. [lfjl. [20l|. There are also mixed 
magneto-optical approaches |6|. For ions, controlled time 
dependent voltages have been used in linear-trap based 
frequency standards [§] and more recently in quantum in- 
formation applications using multisegmented Paul traps 
[2l| , an array of Penning traps [22j , and also in 2D con- 
figurations |23j . 

As said before, an obvious solution, at least in princi- 
ple, to avoid spilling or heating of the atoms is to per- 
form a sufficiently slow (adiabatic) transport. For some 
applications, however, this takes too long. In particular, 
since transport could occupy most of the operation time 
of realistic quantum information alg orithms, "shuttling 
times" need to be minimized @, l2l| . In addition, long 
times may be counterproductive in practice and induce 
overheating from coils or fluctuating fields and decoher- 
ence. There are in summary good reasons to reduce the 
transport time, and indeed several theory and experimen- 
tal works have studied ways to make fast transport also 
faithful [UBS]. 

Invariant-based inverse engineering is ideally suited to 
this end. The main aim of this paper is to set the ba- 
sic invariant-based inverse engineering transport theory, 
analog to the one developed recently for trap expansions 
[27| . We shall also show that previous approaches for 
efficient transport [U, [2f| and some generalizations are 
embraced by it, and point out the potential limitations 
of the method. In Sec. |TT] we shall provide the main 
concepts and formulae of the time dependent quadratic- 
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in-momentum invariants relevant for transport problems. 
The two main reference cases are (z) rigid harmonic os- 
cillator transport and (ii) transport on an arbitrary trap 
with force compensation. In Sec. IIIII we explain and ap- 
ply the inversion technique; this is compared in Sec. IIVI 
with an alternative "bang-bang" approach based on time 
segments of constant acceleration. Section IVl deals with 
practical limitations and Sec. [VI] discusses the results 
and draws the conclusions. 



II. DYNAMICAL INVARIANTS 

In a seminal paper Lewis and Riesenfeld derived a 
simple relation between the solutions of the Schrodinger 
equation of a system with time-dependent Hamiltonian 
and the eigenstates of the corresponding invariants [28| . 
They paid special attention to the time-dependent har- 
monic oscillator and its invariants quadratic in position 
and momentum, related to earlier work by Ermakov on 
the classical oscillator [29| . From a classical physics point 
of view, Lewis and Leach found the general form of the 
Hamiltonian compatible with invariants quadratic in mo- 
mentum [30(, including non harmonic potentials. This is 
the result that interests us here, together with the cor- 
responding quantum formulation by Dhara and Lawande 
31]. In this Section we shall state the main concepts and 
equations and apply them to standard transport prob- 
lems. 

A ID Hamiltonian with an invariant which is quadratic 
in momentum must have the form H = p 2 /2m + V(q, t), 1 
with 0, El 



V(q,t) = -F(t)q + 



m 9 , . , 1 



u 



a(t) 



P(t) 



(1) 



p, a, u>, and F are arbitrary functions of time that satisfy 
the auxiliary equations 



p + u 2 (t)p = ^, 
a + uj 2 (t)a = F(t)/m, 



(2) 
(3) 



with ujq constant. Their physical interpretation in the 
context of transport is detailed below. The dynamical 
invariant, up to a constant factor, is given by 



I = -^—[p(p-ma)-mp(q-a)] 2 
Zm 



U 



q — a 



and verifies 



(4) 



(5) 



1 Following the usual practice, q and p will denote operators or 
numbers, and the context should clarify their meaning. 



so that 4r{'ip(t)\I(t)\tp(t)) = for any wave function ip(t) 
that evolves with H. ip(t) may be expanded in terms of 
constant coefficients c n and eigenvectors ip n of /, 



V'fo*) = ^c„e la '*-0n(<7,£) 

n 

I(t)ip n (q,t) = X n tp n (q,t), 



(6) 
(7) 



where the A n are time-independent eigenvalues. We shall 
generally deal with iji n normalized to one, but continuum, 
delta-normalized states are also possible. The phases a n 
satisfy [H, |H| 



da n 



0_ 

dt 



H 



Ipn ih 
1 

^ ~hj - \p 2 ' 2p 2 
The ip n are in practice obtained easily as [311 ] 

q - 



p 



(8) 
(9) 

(10) 



from the solutions ^(cr) (normalized in er-space) of the 
auxiliary stationary Schrodinger equation 



h 2 d 2 

2m da 2 



1 2 2 

2 



U(a) 



4>n = A„0„. (11) 



Whereas trap expansions and contractions imply a time 
dependent p function [27|, a large family of transport 
problems may be described by taking 



P(t) = l, U 2 (t)=LU 2 



(12) 



so the auxiliary Eq. ([2]) plays no role and only Eq. ([3]) is 
relevant. Except in the final discussion we shall assume 
that the conditions (|12[) hold from now on, and consider 
in detail two main reference cases. 



A. Main cases 

(i) Rigid harmonic oscillator driven by the "transport 
function" qo(t). (Hereafter "harmonic transport" for 
short.) Suppose that a harmonic trap is moved from 
qo(0) at time t = to d = qo{tf) at a time i/. In Eq. ((T|) 
this case corresponds to 



mw$q (t), uj(t) = w , U = 0. 



(13) 



Adding to V the irrelevant time dependent global term 
towq(7q/2, which produces no force, the trap potential 
can be written as a moving harmonic oscillator muj 2 [q — 
qo(t)?/2, 



H = p 2 /2m + mul[q - q (t)} 2 /2, 



(14) 
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and a may be identified with a classical trajectory q c 
since Eq. ^ becomes 

q c + Uo{q c - go) = 0. (15) 

The invariants and transport modes will depend on it. 
In this case A n = E n = (n+ 1/2)?vu)q, and the transport 
mode e lan ip n takes a physically transparent form, 

e M "Vn = e L Jo 2 \e m ^ h Mq-Qc). (16) 

Efficient transport will be engineered in the following sec- 
tion by designing first an appropriate classical trajectory 
q c {t), from which the trap motion trajectory qo(t) is de- 
duced via Eq. (TT5|) . 

A variant of this case is vertical transport with a grav- 
ity force, so that F — muj^qo — mg and Eq. (Tl5j) is 
modified to 

q c + Uo(q c -q ) = -g. (17) 

(») Arbitrary-trap driven transport with compensat- 
ing force. (Hereafter "compensating force approach" for 
short.) Now, in Eq. ^ 

uj = uj = 0, (18) 
F = mq . (19) 

In this case the trap potential U[q — qo(t)] is arbitrary 
(in particular it could be harmonic) , and it is rigidly dis- 
placed along qo(t), so a in Eq. ([3]) may be now identified 
with the transport function go- In addition to U, there 
is a time dependent linear potential term —mqqo in H, 

H = p 2 /2m - mqqo + U(q - q Q ). (20) 

The corresponding force compensates exactly the iner- 
tial force due to the trap motion in the rest frame of the 
trap, so that the wave function in that frame is not mod- 
ified up to a time dependent global phase factor, see the 
Appendix [21 This Hamiltonian has been proposed by 
Mashuda and Nakamura following a very different route, 
using a "fast- forward" scaling technique |25| . 



III. INVERSE ENGINEERING METHOD 

The Lewis-Riesenfeld theory of invariants has been 
considered before in harmonic-oscillator driven transport 
in the direct (rather than inverse) way, by setting the 
transport function go and analyzing the final heating, in 
particular in adiabatic or quasi- adiabatic regimes |2J. We 
shall use it instead as the basis for an inverse engineering 
approach, including also non-harmonic driving. The two 
main cases discussed above require different implementa- 
tions. In both cases we shall assume that go is displaced 
from to d in a time t / . 



A. Harmonic transport 

In case (i) we may adopt as in (27| an inverse engi- 
neering strategy by designing first the classical trajectory 
q c to assure that the transport modes coincide with the 
eigenvectors of the instantaneous Hamiltonian at initial 
and final times. This amounts to impose the commuta- 
tivity of I(t) and H(t) at t — and t — tf, which can be 
achieved by setting, see Eq. (JTHJ) , the following boundary 
conditions, 

?c(0) = <Zo(0) = 0; g c (0) = 0; ? c (0) = 0, (21) 
Qc(t f ) = q (t f ) = d; q c (t f ) = 0; q c (t f ) = 0, (22) 

the last conditions in each line being determined by con- 
sistency with Eq. (fT5l) . q c (t) is then interpolated by 
assuming some flexible functional form, such as a poly- 
nomial, q c (t) = X)n=o^"^™' wnere the /3 n are found by 
solving the system of equations established by the bound- 
ary conditions. The resulting q c depends on time only 
trough the ratio s = t/tf, and is directly proportional to 
d, 

g c (f)/d = 10s 3 -15s 4 + 6s 5 . (23) 
Once q c is fixed wc get the trap trajectory from Eq. (|15p . 

«fo(t)=&(t)/w§ + ffc(t). (24) 

This procedure is equivalent to the one followed by Mur- 
phy et al. [24]], who used a Fourier sum as interpolating 
function for q c . 

For short times tf, see Fig. 1, the corresponding trap 
trajectories qo(t) could exceed the interval [0, d]. For the 
polynomial ansatz (|23l) it occurs symmetrically at lower 
and upper edges of the interval for tf < 2.505/wo « 0.4 
To, where To = 27r/wo is the oscillation period. This may 
or may not be a problem depending on the geometrical 
constraints of the experimental setting. 

An interesting generalization is to consider boundary 
conditions for stopping atoms when their initial average 
velocity is known. Suppose that an atom gun or pulsed 
valve sends atoms with a specific average velocity vq, as in 
coil-gun experiments with paramagnetic atoms [32, [33[ , 
or a Stark decelerator for polar molecules |34|-36j. The 
valve opening time is controlled, so a traveling harmonic 
trap wrapping the atoms can be turned on at time t — 
and moved along some trajectory qo(t) to stop them 
at a fixed distance d in a specified time tf. The final 
conditions may still be given by Eq. (|2"2"]) . but a different 
set of initial conditions are to be imposed, 

q c (0) = 0, j c (0) = v , q c (0) = 0. (25) 

In this case the n-th initial transport mode at t = 
does not coincide with the n-th stationary eigenstate of 
H(0) but with its moving version, e lmv ° q / h 4> n (q). ujq is in 
principle arbitrary, but it may be optimized taking into 
account the spatial width of the incoming state. Its value 
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FIG. 3: (Color online) Stopping atoms: The shaded areas 
represent the values of a = votf/d and b/2n = tf/To for 
which there is a t in (0, tf) where qo(t) < (left figure) or 
qo(t) > d (right figure). We have used the polynomial ansatz 
ofEqs. (125157)1 . 



FIG. 1: (Color online) qo/d versus s — t/tf for tf = 
12.57/wo = 2T (blue line with dots), tf = 2.505/o;o (this 
is the critical value; red line with triangles), and tf = 2/wo 
(brown line with squares). For all three cases q c /d is the 
dashed line, hardly distinguishable from qo/d for the slowest 
case, i.e., for tf = 12.57/cuo- 
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FIG. 2: (Color online) Three stopping trajectories, a = 0.8. 
b = Wot/ = 1.6, tf/To = 0.25 (brown line with squares); 
6 = cjotf = 1.9, tf/To ~ 0.30 (red line with triangles); b — 
uot f = 3.0, t f /To = 0.477 (blue line with circles). They all 
share the same q c (s)/d (black dashed line) given by Eq. ((26j . 



also has an impact on the domain of the trajectory as we 
shall see. The polynomial ansatz gives now 



q c (t) = d[3(2 ~ a)s 5 - (15 - 8a)s 4 
- 2(-5 + 3a)s 3 + as], 

q (t) = d\3{2-a)s 5 - (15 - 8a)s 4 



- 2(-5 + 3a) 



60(2 -a) 



-12 



b 2 

(-5 + 3a) 
b 2 



(26) 



12 (15-8«) s2 



b 2 



(27) 



where s = t/tf, a = vgtf/d and b = uiotf = 2irtf/To. 
Some trajectories are shown in Fig. [2j The shaded re- 
gions in Fig. [3] correspond to the values of a and 6/27T 



for which the trajectory qo(t) exceeds the domain [0,d]. 
Even though the details are now more complicated than 
for the rest-to-rest case, two simple general rules can be 
drawn: qo is never negative when tf > To (the asymptotic 
threshold for large a is at b = 6), whereas if a > 2.513 
there is always some t in (0, tf) for which qo(t) > d. 

Launching or "catapulting" atoms at rest to end up 
with a chosen velocity t>o [lfj may be designed similarly 
by setting Eq. ([2Tj) for the initial conditions, and the 
final boundary conditions as 



q c (t f ) = d, q c (t f ) = v , q c {t f ) = 0. 



(28) 



A major practical concern in all these applications should 
be to keep the harmonic approximation valid. This may 
require an analysis of the actual potential in each specific 
case and of the excitations taking place along the non- 
adiabatic transport process. Without performing such 
detailed analysis, the feasibility of the approach for a 
given transport objective set by the pair d, tf can be 
estimated rather simply by comparing lower excitation 
bounds provided in Sec. fVl with the trap depth. 



B. Compensating force approach 

In the compensating- force approach, case (it), we may 
proceed similarly but now the variable a of Eqs. ([TJ and 
([3]) is directly go so, instead of fixing values of q c and 
its derivatives at the boundary times, we shall fix values 
of qo and its derivatives. For simplicity we shall con- 
sider only a rest-to-rest scenario, but the generalizations 
are straightforward. The compensating potential —mqoq 
in the Hamiltonian (f2"U| should vanish before and after 
transport, since the trap remains at rest and q'o = for 
t < and t > tf. To make the ip n coincide, up to a 
global phase factor, with the eigenstates of the Hamilto- 
nian before and after transport, H 
at t = and t = tf, we impose 



p 2 /2m + U(q-q a ), 



«o(0) 
Qo(tf) 



9o (0) = 0, 
d, qo(tf) - 



(29) 



We may as well impose ijo — as a boundary condition 
at t = and it, to have a continuous qo(t) but, at least 
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formally, it is not strictly necessary. In practice several 
experiments have been designed with (approximate) dis- 
continuities in the trap acceleration [lg, [l^]. We shall 
come back to this point below. Clearly the implementa- 
tion of the compensating-force approach is subjected to 
different limitations from the ones applicable to the har- 
monic transport without compensation. The main prob- 
lem now is not anharmonicity, which is included in the 
theory from the start by admitting an arbitrary U, but 
the feasibility of the compensating force term. According 
to the mean value theorem, see Sec. |Vl a lower bound 
for the maximum of the absolute value of qo is Id J A. 



IV. BANG-BANG ACCELERATION METHODS 



V. TRANSIENT ENERGY EXCITATIONS 

Whereas, ideally, non-adiabatic faithful transport can 
be performed for arbitrary transport distances and times, 
in practice the process could be limited, apart from the 
geometrical constraints discussed in the previous section, 
by the maximal transient excitation energies allowed to 
neglect the effect of anharmonicities of the actual poten- 
tial in case (i), or by the difficulties to implement strong 
compensating forces in case (ii). We shall analyze these 
effects from the point of view of different bounds ob- 
tained for the average (potential) energy using the Euler- 
Lagrange equations, and for the instantaneous potential 
energy and acceleration by means of the mean value the- 
orem. 



It should be clear from section IIII Al that there are in- 
finitely many functions qo(t) which, for harmonic trans- 
port, lead to the ideal boundary conditions. A somewhat 
extreme case, which has however been implemented ex- 
perimentally because of its relative simplicity, is to com- 
bine time segments with a constant acceleration [l6l . [l9| . 
The simplest trap trajectory of this type implies a con- 
stant positive acceleration Ad/t'j from to tf/2, and a 
deceleration —Ad/t'j from tf/2 to tf. The resulting qo is 
formed by two parabolas matched at tf/2, 



Qo/d 



2s 2 , 
4(s 



< s < 1/2 
1/2 < s < 1 



(30) 



The corresponding velocity qo increases linearly from 
to 2d/tf and then decreases from there to 0. 

The classical trajectory q c satisfying Eq. (fT5j) with Eq. 
(f3"0")) . and the boundary conditions ([2~Tj) at t = is given 
by 



qo = 



-4d 
/ - 



id 



(1 



COS Wot) 
f- COS UJ(yt 



2 cos 



LJ 



t - -t- 

L 2 



( 31 ) 

for the first and second time segments. From this result 
one can check that the boundary conditions at tf (|2"2")l are 
satisfied only for a discrete set of times tf t x = AttN/ojo, 
with N <E N, i.e., for multiples of two oscillation periods. 
For all other times this scheme will heat the atoms. A 
perturbation theory analysis shows that, even for the se- 
lected discrete times, the bang-bang method is less stable 
than the inverse invariant method with respect to an an- 
harmonic perturbation of the transporting trap potential. 
The details are shown in Appendix [B] 

A variant of this method, using e.g. the trap trajec- 
tory (f3"U|) and, in addition, compensating forces as in case 
(ii), may be appealing in practice because it is relatively 
simple to implement the compensating force, at least ap- 
proximately, as a piecewise function. 



A. Quasi-optimal trajectories for harmonic 
transport 

The instantaneous average energy for a harmonically- 
driven transport mode can be calculated from Eqs. 

mm, 



(^ n (t)\H(t)\M)) 
1 



hwo 



Tfl , o 1 9 / \ 2 

., , ■ y9 c + 2 ma Wc-<7o) ■ 



(32) 



Moreover the instantaneous average potential energy is 
(V(t)) = (n + l/2)+E P . The first, "internal" contri- 
bution remains constant for each n, and the second term 
has the simple form of a potential energy for a classical 
particle, Ep = ^muJo(q c — qo) 2 ■ Its time average, using 
the relation (|15l) between qo and q c , takes the form 



E P = 



2t f wl 



q 2 dt. 



(33) 



We can use a generalized Euler-Lagrange equation 
d^qc/dt 4 = to minimize this integral subject to four 
boundary conditions [37| . the ones for q c and q c in Eqs. 
(|21I22|) . This results in a "quasi-optimal" classical tra- 
jectory 



q c = d(3s 2 - 2s 3 ). 



(34) 



Whereas Eq. (|34|) does not satisfy the six boundary con- 
ditions (|21l22j) . it provides in any case a lower bound for 
the time average of Ep, as the set of functions satisfying 
the six conditions is smaller than the one satisfying four 
of them. Substituting Eq. (fMf into Eq. (|3"3")l one gets 
the desired lower bound, 



E P > 



6md 

+4, ,2 



(35) 



This bound describes the relevant dependences, as shown 
by numerical comparisons with actual time-averaged en- 
ergies for polynomial trajectories, Eq. (|2"3")1 . see Fig. U{a), 
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FIG. 4: (Color online) Dependences of time-averages energies 
with tf. The 5 is the asymptotic exponent oft/. Parameters: 
uiq = 2-7T x 8 Hz, _Eo = huo/2, d = 2.25 mm, and rubidium- 
87 atoms that begin and end in the ground vibrational state, 
(a) Bound ()35|) (solid red line), and time average of Ep for 
a polynomial trajectory (dotted blue line), (b): AA Bound 
(dotted-dashed magenta line), and time average of AH for a 
polynomial trajectory (dotted blue line). 



orthogonal initial and final states [3|| Eo| is AHtf > 
fi,arccos[exp(— mujod 2 /4h)], which tends to the result for 
orthogonal states when d >> (Ah/majo) 1 / 2 .) As it oc- 
curs for harmonic trap expansions [38| . while certainly 
correct as a bound, it does not describe the dependences 
found for the averaged standard deviation for fast pro- 
cesses (small tf). 



B. Mean value theorem 



slope x„ 




and sets a rather strong tf scaling, compare this with 

the milder dependence on t J 2 of the time-averaged tran- 
sient energy in invariant-based inverse-engineered expan- 
sion processes [H| . 

As for the variance (AH) 2 = (ipn\H 2 \'if) n )-(ij n \H\'il) n ) 2 
for the n-th transport mode, it takes, after a somewhat 
lengthy calculation, a simple form, 



(AH) 2 = 2hu (n + l/2) 



moJoilc - qo) 2 



1 



-mq c 



(36) 

Using again an Euler-Lagrange equation we find for 
its time average the lower bound (AH) 2 > 12h(n + 
l/2)md 2 /uJot A f. This does not establish a lower bound 
for the average of the standard deviation AH, but 
agrees with the scaling with tf that we observe nu- 
merically as tf -> 0, AH cx tj 2 , see Fig. HJb). 
This should be contrasted with the Aharonov-Anandan 
(AA) relation [39] AH > h/Atf, applied to transport 
among orthogonal states. (The general expression for 
ground-state to ground-state transport allowing for non- 



FIG. 5: Graphical representation of the lower bound of 
Eq. (f57) for the maximum velocity x ma x (upper graph), and 
maximum acceleration |i| ma i of Eqs. ((38)) and (|39l) (lower 
graph). 

The mean value theorem (MVT) sets another useful 
bound since it applies to the instantaneous values rather 
than to a time average. The argument may be applied to 
q c in case (i) or to go i n case (ii), so we shall formulate 
it in terms of a generic x, assumed continuous in [0,f/], 
differentiable in (0, tf), and such that x(0) = x(0) = 
x(tf) = and x(tf) — d. The maximum of its time 
derivative must be 



> d/t f 



(37) 



at some point t m in (0,tf), see Fig. [5j We can now 
use that point to divide [Q,t/] into two segments and 
apply the MVT again, now to the derivative. In the first 
segment the derivative goes from to x max so 



l^lmax ^ dj (tftjyi). 



(38) 
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Similarly in the second segment, 

\x\ma X >d/[t f (tf-t m )], (39) 

and, for the whole interval, |x| mQ:E > d/[tf min(i m , tf — 
t m )]. Irrespective of the location of the point t m , 
min(t m ,tf —t m ) < tf/2. We can thus set a lower bound 
for the (absolute value of) the maximum of the second 
derivative, 

\x\rnax > U/t). (40) 

When x = q c , case (i), Eq. (|40p gives a bound for the 
instantaneous quasi-potential energy 

E P >2m(^j . (41) 

In this case, however, we get a tighter lower bound di- 
rectly from the time average (1351) . 

With x = qo, case (it), we get a lower bound for the 
maximum trap acceleration. In particular the accelera- 
tions of compensating forces, typically limited by exper- 
imental constraints such as gradient magnetic fields or 
Stark electric fields, should at the very least reach the 
value 2d/t 2 f. 

VI. DISCUSSION 

We have applied the Lewis-Riesenfeld method [28| for 
quadratic-in-p invariants [3ll | , combined with inverse en- 
gineering of trap trajectories, to design fast and faithful 
atomic transport. The limitations have been quantified, 
and relations to other approaches that can be included in 
this framework have been pointed out. Another approach 
to accelerate adiabatic processes has been recently pro- 
posed by Berry [4JJ and can as well be applied formally to 
transport. The central idea is to construct a Hamiltonian 
that drives the system exactly along the adiabatic ap- 
proximation defined for some reference time-dependent 
Hamiltonian Ho(t), without transitions among the in- 
stantaneous eigenstates of Hq. For an arbitrary, rigid 
trap potential moving from to d along a path qo(t), 
the eigenvectors of the instantaneous (reference) Hamil- 
tonian Hq = p 2 /2m + U(q — qo) are simply displaced from 
the original location, 

\n(t)) =e- ipqo{t)/h \n(Q)). (42) 

The transitionless driving Hamiltonian Hq + 
ih^2 n \dn)(n\ becomes in this case Hq + pqo, com- 
pare with Eq. (|2"0"|) . Note the freedom to choose Hq. It 
could even be suppressed during transport: the simple 
Hamiltonian pqQ would also keep the same populations 
(of any Hq) without generating dynamical phases e~ lEnt 
for each eigenvalue. This is a rather intuitive result since 
the corresponding propagator is nothing but the dis- 
placement operator e~ %pqot ^^ h . This approach provides 



thus a formal solution to fast and faithful transport, 
but the practical realization of a qap Hamiltonian term 
remains an open question. 

With respect to the general framework embraced by 
Eqs. HEI]), the studied cases termed (t) and (it) in 
the main text are very relevant but not exhaustive. For 
U = 0, the most general case occurs when transport is ac- 
companied by expansions and contractions, so that the 
time dependence of uj(t) and p(t) has to be considered 
if the invariant-based inverse engineering method is ap- 
plied. 

For harmonic trap expansions or contractions in the 
gravity field (q becomes a vertical coordinate), F — —mg 
and the center of the trap suffers a time dependent trans- 
lation go = —2g/ui 2 . Again, a in Eq. © may be inter- 
preted as a classical trajectory, now of a time dependent 
harmonic oscillator subjected to the gravity field, 

q c + u 2 (t)q c = -g. (43) 

This is admittedly not a proper transport problem, but 
its formal treatment is the same, and has recently been 
implemented experimentally [42| . also for Bose-Einstein 
condensates [43j]. 

As for further extensions or open questions of the in- 
variant approach, one may investigate the use of more 
complex invariants, not restricted to being quadratic in 
p [4J], in particular to tackle anharmonic transport. In 
the frame of quadratic-in-p invariants, anharmonic traps 
can be dealt with by a compensating force (case (it)), 
but this force might be difficult to implement for large 
accelerations. If F = 0, uj = 0, and U(q — a) ^ in 
Eq. (JTJ) , a should be the trap trajectory q Q , which is only 
consistent with the auxiliary equation a — 0, see Eq. ([3]), 
for constant-velocity trajectories, incompatible with the 
boundary conditions (|2"9"|) . A way out, to be explored, 
may be to use the invariants to implement minimization 
algorithms of the final vibrational excitation. 

Finally, further work will be devoted to understand 
and mitigate the effects of noise, for which the inversion 
method is intrinsically robust [24| , and of atom-atom in- 
teractions. The regime of Tonks- Girardeau gas can be 
treated along similar lines of the single particle case and 
for condensates, scaling techniques may be used as in 
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Appendix A: Displacement unitary transformation 



quartic term as a perturbation, 
H = H + V u 

[q-qo(t)} 2 



2m 



x 



R 



Vi = -V a 



[q-Qo(t)Y 



V 



[Q-Qo(t)] 4 



(B3) 



Consider the following (time dependent) position and 
momentum displacement unitary operator 



U 



e ipq {t) / h e -imq {t)q/ H 



(Al) 



We could as well use variants of U with different order- 
ings without affecting the final result. Starting from the 
Schrodinger equation 



mtW) = H\ij)), 

where, as in Eq. P0]l. H = p 2 /2m + U{q - q a ) 
the corresponding equation for |4>) = li\tp) is 

ihd t \<f>) = UHlfi\<$>) +ih(d t U)U f \<l>) 



2m 



U(q) 



mq 21 



(A2) 
mqq , 

(A3) 



Any stationary state in this "trap frame" will remain so 
in spite of the trap motion thanks to the compensating 
effect of the term —mqqo in H. 



Appendix B: Perturbation theory analysis of the 
effect of anharmonicity 

In this Appendix we shall use perturbation theory to 
determine the effect of small anharmonicities using in- 
verse or bang-bang trap trajectories. We start for con- 
creteness from the "cigar trap" potential associated with 
a Gaussian beam with a moving focus, 



V(A,r) = -V e- 2r2 / w2 ^—!- K 

1 + 



(Bl) 



where r and A are radial and longitudinal coordinates, 
A = q-q (t), 



;(A) =io Wl + 



A_ 



(B2) 



is the spot size, xr — nw 2 /\ is the Rayleigh length, and 
wq the waist. 

For a tight radial confinement we may ignore the radial 
coordinate and set r = 0. The resulting longitudinal po- 
tential can be expanded around the minimum. Retaining 
the first correction to the harmonic term, and ignoring 
the constant, we split the Hamiltonian considering the 



where Vq = mLj^x 2 R j1. Using time dependent perturba- 
tion theory we calculate the overlap (V>(if )|V>(i/)) be- 
tween the state evolving with the harmonic oscillator 
\tp(t)) and the perturbed state \ip(t)) at the final time 
t = tf. \i/>(t)) is chosen as the transport mode (|16l) . with 
(I2T1) satisfied. 

We approximate the perturbed state in first order as 

\m) = u (t, o)i^(o)) 

- l -J^ dt'U Q {t,t')V{t')U {t\Q)${Q)),{m) 



where 



so 



£/ (t,0) =exp 



^ Jo 



dt'H (t') 



(B5) 



m f )m f )) = m f )\u (t f ,Q)\i>(o)) 



(B6) 



dt'(i;(t f )\Uo(t f ,t')V(t')Uo(t',0M(0)). 



At t = the initial state is also an eigenstate of the 
harmonic oscillator, \ip(0)) = |^(0)), so the first term in 
the right side of the previous equation is 1. Using the 
transport modes in Eq. (|16p . we calculate the bracket 
term in Eq. (|B6[) . 



(B7) 



m f )\uo(t f ,t')v(t')Uo(t',o)\m} 



-V a 



m'^q-q^t'tm')) 



Using for q$ and q c the functions for the bang-bang case 
described in Sec. |lVl see Eqs. ((30)) and (|3T|) . and per- 
forming the time integral we arrive, in first order, at 



Wf)Mt f )) = i--r ) 



(B8) 



where the result for T is explicit but lengthy and not very 
illuminating so it is omitted here. It may be checked that 
when tf — > then T — > 0, but this is not a very useful 
limit since the bang-bang procedure will not work for 
times smaller than Ait /luq . As a more useful case we take 
for tf the discrete times tf^ — AttN/ojo with TV G N [l9l |. 
Then J- takes the form 



_2-(io+n)( 2n )!] 
N 7 mir 7 n\u)QX 2 R 



(B9) 



x {l5367V 8 ft 2 [l + 2n{l + n)]w 8 

+ 576d 2 N 4 hm(l + 2n)TT 4 u} + 35d' l m 2 u} 2 }. 
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In the limit xr — > oo, Tbb — > 0. Increasing the waist, and 
keeping the other parameters constant, the potential is 
more harmonic. By contrast, as xr — > 0, J-bb — oo. 
Also, !Fbb —oo when luq — > and luq — > oo. 

If instead of the bang-bang functions we choose the 
inverse-engineered polynomial q c in Eq. (|23p and the 
corresponding qo, T becomes, in the general case, i.e., 
for arbitrary parameters and in particular an arbitrary 

*/> 



J~ inn — 



_ 2 -(™+ 2 )Vo(2n)!! 



1728000d 4 
lOOli^w 8 , 



1440d 2 S,(l + 2n) fr 2 [3 + 6n(l + n)t f ] 



For the final times t/.^v, 



. 2 -(8+»)3(2n)!! 



(BIO) 



1001N 7 mir 7 n\v x R 
x {l281287V 8 /i 2 [l + 2n(l + n)]7r 8 

+ 34320d 2 iV 4 /im(l + 2n)n 4 uj + 1125d 4 m 2 w 2 } 



Comparing the factors in (|B9[) and (|B10|) . we see that 
\Finv\ < \Fbb\ for > 0. 
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